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Introduction 


This report summarizes work accomplished under NASA Grant NCC2-5 167, 
"Development of Advanced Methods of Structural and Trajectory Analysis for Transport 
Aircraft," during the first year of the grant October 1, 1995 - September 30, 1996. The 
effort was in two areas: (1) development of advanced methods of structural weight 
estimation, and (2) development of advanced methods of flight path optimization. 

During the Spring of 1996 both graduate student research assistants working on the 
project, H.C. Chou and Mark Chambers, resigned to take positions in industry. This 
required assigning three new Santa Clara people to the project: Dr. Lee Homberger, 
Associate Professor of Mechanical Engineering; Robert Windhorst, graduate student 
research assistant; and Frank Dickerson, undergraduate student. These new people 
inevitably required time to learn the ACSYNT code and the nature of the ongoing 
research. The result is that some of the tasks in the work statement have not been 
completed at this time, but will be at the completion of the Grant. 

Dr. H. Miura and M. Moore were the NASA collaborators on the Grant. 


Review of Results in Structural Weight Estimation 

A report that was prepared under a previous grant was published in May 1996 
("Analytical Fuselage and Wing Weight Estimation of Transport Aircraft," by 
M. Ardema, M. Chambers, A. Patron, A. Hahn, H. Miura, and M. Moore, NASA TM 
1 10392). A paper that summarizes this report has been accepted for presentation at the 
World Aviation Congress to be held in October 1996; a copy of this paper appears in 
Appendix A. 

Throughout the year, integration of the structural weight computer code, PDCYL, into 
ACSYNT has continued. Input variables used by PDCYL but already in ASCYNT have 
been removed from PDCYL. Infrequently used input variables have been defaulted. Data 
transfer has been modified so that optimization runs with ACSYNT can be done with 
PDCYL as an integral part of the code. 

The major effort to the first year of the grant was to develop an improved method of 
estimating the weight of wing and fuselage structures made from composite materials. 

This involved an extensive literature search, the coding of a composite materials 
subroutine, and demonstrating the code. This work is discussed in detail in Appendix B. 

Previously in ACSYNT, the weight of composite material structures was estimated 
assuming quasi-isotropic materials, maximum stress failure theory, and smeared structural 
elements. The capability being developed accounts for realistic lay-ups of unidirectional 
fiber/matrix composites and uses a bi-axial strain failure theory. The new composite 


routine has been implemented for the fuselage weight calculation and will be implemented 
for the wing weight calculation in the second year of the grant. 

A user’s manual for the new composite subroutine may be found in Appendix B. As a 
check case for the new subroutine, the weight of a composite fuselage of the ASA 2150 
has been estimated. 

The final effort in the structures area has been support of the project to design and 
analyze a 150 passenger advanced transport airplane, the ASA 2150. PDCYL has been 
used as an integral part of ACSYNT to estimate the fuselage and wing weights of this 
aircraft. Appendix C gives the details of the weight calculations for both Aluminum and 
Graphite/Epoxy fuselage versions of the ASA 2150. 

Appendix C shows that at a gross take-off weight of 152,181 pounds, the ASA 2150 is 
estimated to have a wing weight of 10,315 pounds and a fuselage weight of 15,652 
pounds when made of Aluminum. Figure 1 shows the ASA 2150 fuselage bending 
moment distribution. The critical loading condition for most of the fuselage is either the 
landing condition (L) or the runway bump condition (B), with a small portion governed by 
the maneuver condition (M). The shell unit weight distribution is shown on Figure 2. 
Approximately the first half of the fuselage is sized by minimum gage, with most of the 
rest yield strength critical. 

When the fuselage is made of composite material, the weight is estimated to be 15,375 
pounds, a weight savings of about 2 % relative to aluminum. The composite material is a 
uni-directional tape made from Hercules AS4 carbon fiber in Fiberite 12K/938 resin. The 
reason for this relatively low weight savings is that for relatively small and lightly loaded 
aircraft such as the ASA 2150, the fact that the composite material thickness must be in 
integer thicknesses of the basic stack thickness means that the structure is in many places 
considerably overdesigned. The basic stack used was a quasi-isotropic lay-up of eight 
unidirectional plys. Also, the nonoptimun factor used for the composite was 17% higher 
than that for the Aluminum design As for the Aluminum design, the composite fuselage 
was sized by minimum gage and yield strength. 


Review of Results in Trajectory Optimization 

Because of the unexpected loss of a senior, experienced research assistant, the analysis 
of the altitude jumps in energy climb paths could not be completed. Rather, a new, basic 
look was taken at energy-state and related approximations. This analysis is related in 
Appendix D and will result in an important new addition to the trajectory optimization in 
ACSYNT. 

Previously, the approach was to minimize a weighted sum of time and fuel 
consumption: 


J = £(*,+*, CT)* 


Under energy-state approximation this may be written 


where 


P = 


v(t - D) 


w 


Thus the optimal climb path is given by 


maxf p 
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This analysis, however, is for range not specified. When range is specified, the optimal 
path is determined by 
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Where X x is the adjoint variable associated with the range equation. It is determined 
from 


K x +K 2 C c T c 

~ v 
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Where C c , T c , and V c are associated with the optimal cruise point, obtained by 
maximizing the Brequet factor within the flight envelope 



maxT V(L/D) 
b,E l C . 


A search for the optimal cruise point, and the evaluation of is being added to the 
ACSYNT code. 
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ABSTRACT 

A method of estimating the load-bearing fuselage weight 
and wing weight of transport aircraft based on fundamental 
structural principles has been developed. This method of 
weight estimation represents a compromise between the rapid 
assessment of component weight using empirical methods 
based on actual weights of existing aircraft, and detailed, but 
time-consuming, analysis using the finite element method. The 
method was applied to eight existing subsonic transports for 
validation and correlation. Integration of the resulting com- 
puter program, PDCYL, has been made into the weights- 
cal culaiing module of the AirCraft S YNThesis (ACSYNT) 
computer program. ACSYNT has traditionally used only 
empirical weight estimation methods; PDCYL adds to 
ACSYNT a rapid, accurate means of assessing the fuselage 
and wing weights of unconventional aircraft. PDCYL also 
allows flexibility in the choice of structural concept, as well 
as a direct means of determining the impact of advanced 
materials on structural weight. 

INTRODUCTION 

A methodology based on fundamental structural 
principles has been developed to estimate the load-carrying 
weight of the fuselage and basic box weight of the wing 
for aircraft, and has been incorporated into the AirCraft 
S YNThesis program (ACSYNT). This weight routine is also 
available to run independently of ACSYNT, and is a modifi- 
cation of a collection of previously developed structural 
programs. 1 *" 4 The main subroutine called by ACSYNT is 
PDCYL. This study has concentrated on modem transport 
aircraft because of the detailed weight information available, 
allowing the weights output from PDCYL to be compared to 
actual structural weights. The detailed weight statements also 
allow nonoptimum factors to be computed which, when 
multiplied by the load-bearing structural weights calculated by 
PDCYL, will give good representative total structure weight 
estimates. These nonoptimum factors will be computed 
through a regression analysis of a group of eight transport 
aircraft. 


PDCYL is able to model both skin-stringer-frame and 
composite sandwich shell fuselage and wing box construc- 
tions. Numerous modifications werfif made to PDCYL and its 
associated collection of subroutines. These modifications 
include the addition of detailed fuselage shell geometry 
calculations; optional integration of a cylindrical tuselage 
midsection between the nose and tail sections; addition of 
landing and bump maneuvers to the load cases sizing the 
fuselage; ability to introduce an elliptical spanwise lift load 
distribution on the wing; variation of wing thickness ratio from 
tip to root; ability to place landing gear on the wing to relieve 
spanwise bending loads; distribution of propulsion system 
components between wing and fuselage; and the determination 
of maximum wingtip deflection. 

BRIEF DESCRIPTION OF ACSYNT 

The Aircraft Synthesis Computer program, ACSYNT, is 
an integrated design tool used in the modeling of advanced 
aircraft for conceptual design studies ACSYNT development 
began at NASA Ames Research Center in the 1970s and 
continues to this day. The ACSYNT program is quite flexible 
and can model a wide range of aircraft configurations and 
sizes, from remotely piloted high altitude craft to the largest 
transport 

The ACSYNT program uses the following modules, not 
necessarily in this order Geometry, Trajectory, Aerody- 
namics, Propulsion, Stability, Weights, Cost, Advanced 
Aerodynamic Methods, and Takeoff. An ACSYNT run would 
normally progress as follows: the Geometry module is called 
to define the aircraft shape and configuration; the Trajectory 
module then runs the vehicle through a specified mission, 
finally the Weight and Cost modules are executed. To deter- 
mine the performance of the vehicle at each mission poinu the 
Trajectory module will call the Aerodynamics and Propulsion 
modules. 


* Work of the first two authors was supported by NASA Ames 
Research Center Grant NCC2-5068. 
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After the mission is completed, the calculated weight of the 
aircraft may be compared with the initial estimate and an 
iteration scheme run to converge upon the required aircraft 
weight This process is necessarily iterative as the aircraft 
weight ACSYNT calculates is dependent upon the initial 
weight estimate. 

ACSYNT is able to perform a sensitivity analysis on any 
design variable, such as aspect ratio, thickness-to-chord ratio, 
fuselage length or maximum fuselage diameter Sensitivity is 
defined as (change in objective funciion/valuc of objective 
function) divided by (change in design variableAJesign 
variable). As an example, if gross weight is the objective 
function and decreases when the wing thickness-to-chord ratio 
increases, then the sensitivity of Lhickness-to-chord ratio is 
negative. It is important to note that while this increase in 
thickness-to-chord ratio lowers the gross weight of the aircrafts 
it may also have a detrimental effect on aircraft performance. 

ACSYNT is also able to sue multiple design variables by 
optimizing the objective function The objective function 
represents the interactions between design disciplines such as 
structures, aerodynamics and propulsion. The automated 
sizing of design variables during the optimization process is 
accomplished using the gradient method Two types of 
constraints may be imposed during the optimization process. 
These are performance -based constraints such as runway 
length or maximum roll angle, and side constraints on design 
variables such as limitations on wing span or fuselage length. 
ACSYNT never violates uvistraints during the optimization 
process so that each iteration produces a valid aircraft. 

METHODS OF WEIGHT ESTIMATION 

Two methods are commonly available to estimate the 
load-bearing fuselage weight and wing box structure weight of 
aircraft. These methods, in increasing order of complexity and 
accuracy, are empirical rcgrevMtu and detailed finite element 
structural analysis. Each method has particular advantages and 
limitations which will be botfly discussed in the following 
sections. There is an additxnal method based on classical 
plate theory (CPT) which may be used to estimate the weight 
of the wing box structure. 

EMPIRICAL 

The empirical approach is the simplest weight estimation 
tool. It requires knowledge of fuselage and wing weights from 
a number of similar existing aircraft in addition to various key 
configuration parameters of these aircraft in order to produce a 
linear regression. This regression is a function of the configu- 
ration parameters of the existing aircraft and is then scaled to 
give an estimate of fuselage and wing weights for an aircraft 
under investigation. Obviously, the accuracy of this method is 
dependent upon the quality and quantity of data available for 
existing aircraft. Also, the accuracy of the estimation will 
depend on how closely the existing aircraft match the configu- 
ration and weight of the aircraft under investigation. All of the 
empirical regression functions currently in the ACSYNT 
program give total fuselage weight and total wing weight. 


FINITE ELEMENT 

Finite element analysis is' the matrix method of solution of 
a discretized model of a structure. This structure, such as an 
aircraft fuselage or wing, is modeled as a system of elements 
connected to adjacent elements at nodal points. An element is 
a discrete (or finite) structure that has a certain geometric 
makeup and set of physical characteristics. A nodal force acts 
at each nodal point, which is capable of displacement. A set of 
mathematical equations may be written for each element 
relating its nodal displacements to the corresponding nodal 
forces. For skeletal structures, such as those composed of rods 
or beams, the determination of element sizing and correspond- 
ing nodal positioning is relatively straightforward. Placement 
of nodal points on these simple structures would naturally fall 
on positions of concentrated external force application or 
joints, where discontinuities in local displacement occur. 

Continuum structures, such as an aircraft fuselage or 
wing, which would use some combination of solid, fiat plate, 
or shell elements, are not as easily discretizable. An approxi- 
mate mesh of dements must be made to model these struc- 
tures. In effect, an idealized model of the structure is made, 
where the element selection and sizing is tailored to local 
loading and stress conditions. 

The assembly of elements representing the entire structure 
is a large set of simultaneous equations that, when combined 
with the loading condition and physical constraints, can be 
solved to find the unknown nodal forces and displacements. 

The nodal forces and displacements are then substituted back 
into the each element to produce stress and strain distributions 
for the entire structural model. 

CLASSICAL PLATE THEORY 

CPT has been applied to wing structure design and weight 
estimation for the past 20 years. Using CPT a mathematical 
model of the wing based on an equivalent plate representation 
is combined with global Ritz analysis techniques to study the 
structural response of the wing. An equivalent plate model 
does not require detailed structural design data as required for 
finite element analysis model generation and has been shown 
to be a reliable model for low aspect ratio fighter wings. 
Generally, CPT will overestimate the stiffness of more 
flexible, higher aspect ratio wings, such as those employed 
on modem transport aircraft. Recently, transverse shear 
deformation has been included in equivalent plate models to 
account for this added flexibility. This new technique has been 
shown to give closer representations of tip deflection and 
natural frequencies of higher aspect ratio wings, although it 
still overestimates the wing stiffness. No fuselage weigh! 
estimation technique which corresponds to the equivalent plate 
model for wing structures is available. 

NEED FOR BETTER, INTERMEDIATE METHOD 

Preliminary weight estimates of aircraft are traditionally 
made using empirical methods based on the weights of 
existing aircraft, as has been described. These methods, 
however, are undesirable for studies of unconventional aircraft 
concepts for two reasons. First since the weight estimating 
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formulas arc based on existing aircraft, tbeir application to 
unconventional configurations (i.e., canard aircraft or area 
ruled bodies) is suspect. Second, they provide no straight- 
forward method to assess the impact of advanced technologies 
and materials (i.e., bonded construction and advanced 
composite laminates). 

On the other hand, finite-element based methods of 
structural analysis, commonly used in aircraft detailed design, 
are not appropriate for conceptual design, as the idealized 
structural model must be built off-line. The solution of even a 
moderately complex model is also computationally intensive 
and will become a bottleneck in the vehicle synthesis. Two 
approaches which may simplify finite-element structural 
analysis also have drawbacks. The first approach is to create 
detailed analyses at a few critical locations on the fuselage and 
wing, then extrapolate the results to the entire aircraft, but this 
can be misleading because of the great variety of structural, 
load, and geometric characteristics in a typical design. The 
second method is to create an extremely coarse model of the 
aircraft, but this scheme may miss key loading and stress 
concentrations in addition to suffering from the problems 
associated with a number of detailed analyses 

The fuselage and wing structural weight estimation 
method employed in PDCYL is based on another approach, 
beam theory structural analysis. This results in a weight 
estimate that is directly driven by material properties, load 
conditions, and vehicle size and shape, and is not confined to 
an existing data base. Since the analysis is done station-by- 
station along the vehicle longitudinal axis, and along the wing 
structural chord, the distribution of loads and vehicle geometry 
is accounted for, giving an integrated weight that accounts for 
local conditions. An analysis based solely on fundamental 
principles will give an accurate estimate of structural weight 
only. Weights for fuselage and wing secondary structure, 
including control surfaces and leading and trailing edges, and 
some items from the primary structure, such as doublers, 
cutouts, and fasteners, must be estimated from correlation to 
existing aircraft. 

The equivalent plate representation, which is unable u> 
model the fuselage structure, is not used in PDCYL. 

METHODS 

OVERVIEW 

Since it is necessary in systems analysis studies to be able 
to rapidly evaluate a large number of specific designs, the 
methods employed in PDCYL are based on idealized vehicle 
models and simplified structural analysis. The analyses of the 
fuselage and wing structures are performed in different 
routines within PDCYL, and, as such, will be discussed 
separately. The PDCYL weight analysis program is initiated at 
the point where ACSYNT performs its fuselage weight 
calculation. PDCYL first performs a basic geometrical sizing 
of the aircraft in which the overall dimensions of the aircraft 
are determined and the propulsion system, landing gear, wing, 
and lifting surfaces are placed. 


Fuselage 

The detailed fuselage analysis starts with a calculation of 
vehicle loads on a station-by-station basis. Three types of 
loads are considered — longitudinal acceleration (applicable to 
high-thrust propulsion systems), tank or internal cabin pres- 
sure, and longitudinal bending moment. All of these loads 
occur simultaneously, representing a critical loading condition. 
For longitudinal acceleration, longitudinal stress resultants 
caused by acceleration are computed as a function of longi- 
tudinal fuselage station; these stress resultants are compressive 
ahead of the propulsion system and tensile behind the propul- 
sion system. For internal pressure loads, the longitudinal 
distribution of longitudinal and circumferential (hoop) stress 
resultants is computed for a given shell gage pressure 
(generally 12 psig). There is an option lo either use the 
pressure loads to reduce the compressive loads from other 
sources or not to do this; in either case, the pressure loads are 
added to the other tensile loads. 

The following is a summary of the methods used; the 
details may be found in Ref. 6. 

Longitudinal bending moment tUstributions from three 
load cases are examined for the fuselage. Loads on the 
fuselage are computed for a quasi-static pull-up maneuver, a 
landing maneuver, and travel over runway bumps. These three 
load cases occur at user-specified tractions of gross takeoff 
weight. Aerodynamic loads are computed as a constant 
fraction of fuselage planform area and are considered negli- 
gible for subsonic transports. For pitch control there is an 
option to use either elevators mounted on the horizontal tail 
(the conventional configuration) or elevons mounted on the 
trailing edges of the wing. The envelope of maximum bending 
moments is computed for all three load cases and is then used 
to determine the net stress resultants at each fuselage station. 

After the net stress resultants are determined at each 
fuselage station, a search is conducted at each station to 
determine the amount of structural material required to 
preclude failure in the most critical condition at the most 
critical point on the shell circumference. This critical point is 
assumed to be the outermost fiber at each station. Failure 
modes considered are tensile yield, compressive yield, local 
buckling, and gross buckling of the entire structure. A 
minimum gage restriction is also imposed as a final criterion. 

It is assumed that the material near the neutral fiber of the 
fuselage (with respect to longitudinal bending loads) is 
sufficient to resist the shear and torsion loads transmitted 
through the fuselage. For the shear loads this is a good 
approximation as the fibers farthest from the neutral axis will 
carry no shear. Also, for beams with large fineness ratios 
(fuselage length/maximum diameter) bending becomes the 
predominant failure mode. 

The maximum stress failure theory Is used for predicting 
yield failures. Buckling calculations assume stiffened shells 
behave as wide columns and sandwich shells behave as 
cylinders. The frames required for the stiffened shells are sized 
by the Shanley criterion. This criterion is based on the premise 
that, to a first-order approximation, the frames act as elastic 
supports for the wide column. 7 

There are a variety of structural geometries available for 
the fuselage. There Is a simply stiffened shell concept using 
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longitudinal frames. There are three concepts with Z -stiffened 
shells and longitudinal frames; one with structural material 
proportioned to give minimum weight in buckling, one with 
buckling efficiency compromised to give lighter weight in 
minimum gage, and one a buckling-pressure compromise. 
Similarly, there are three truss-core sandwich designs, two for 
minimal weight in buckling with and without frames, and one 
a buckling-minimum gage compromise. 

It is assumed that the structural materials exhibit elasio- 
plastic behavior. Further, to account for the effects of creep, 
fatigue, stress-corrosion, thermal cycling and thermal stresses, 
options are available to scale the material properties of 
strength and Young’s modulus of elasticity. In the numerical 
results of this study, all materials were considered elastic and 
the full room -temperature material properties were used. 

Composite materials can be modeled with PDCYL by 
assuming them to consist of orthotropic lamina formed into 
quasi-isotropic (two-dimensionally, or planar, isotropic) 
laminates. Each of the lamina is assumed to be composed of 
filaments placed unidirectionally in a matrix material. Such a 
laminate has been found to give very nearly minimum weight 
for typical aircraft structures. 

Wing 

The wing structure is a multi-web box beam designed by 
spanwise bending and shear. The wing -fuselage carry through 
structure, defined by the wing-fuselage intersection, carries the 
spanwise bending, shear, and torsion loads introduced by the 
outboard portion of the wing. 

The load case used for the wing weight analysis is the 
quasi-static pull-up maneuver. The applied loads to the wing 
include the distributed lift and inertia forces, and the point 
loads of landing gear and propulsion, if placed on the wing. 
Fuel may also be stored in the wing, which will relieve 
bending loads during the pull-up maneuver. 

The wing weight analysis proceeds in a similar fashion to 
that of the fuselage. The weight of the structural box is 
determined by calculating the minimum amount of material 
required to satisfy static buckling and strength requirements at 
a series of spanwise stations. The covers of the multi-web box 
are sized by buckling due to local instability and the webs by 
flexure-induced crushing. Required shear material is computed 
independently of buckling material. Aeroelastic effects are not 
accounted for directly, although an approximation of the 
magnitude of the tip deflection during the pull-up maneuver is 
made. For the carry through structure, buckling, shear, and 
torsion material are computed independently and summed. 

As for the fuselage, there are a variety of structural 
geometries available. There are a total of six structural 
concepts, three with uastiffened covers and three with truss- 
stiffened covers. Both cover configurations use webs that are 
either Z-stiffened, unflanged, or trusses. 


GEOMETRY 

Fuselage 

The fuselage is assumed to be composed of a nose 
section, an optional cylindrical midsection, and a tail section. 
The gross density and fineness ratio are defined as 


W B 

P8= 'fj 

(1) 


(2) 


where Wg is the fuselage weight {Wg = gross takeoff weight 
excluding the summed weight of the wing, tails, wing- 
mounted landing gear, wing-mounted propulsion, and fuel if 
stored in the wing), Vg is the total fuselage volume, lg is the 
fuselage length, and D is the maximum fuselage diameter 
The fuselage outline is defined by two power-law bodies of 
revolution placed back-to-back, with an optional cylindrical 
midsection between them (Fig. I). (For the present study, all 
eight transports used for validation of the analysis used the 
optional cylindrical midsection ) 

The horizontal tail is placed according to its quarter chord 
location as a fraction of the fuselage length. 

Propulsion may be either mounted on the fuselage or 
placed on the wing. In the case of fuselage mounted propul- 
sion, the starting and ending positions of the propulsion unit 
are again calculated from their respective fractions of fuselage 
length 

Similarly, the nose landing gear is placed on the fuselage 
as a fraction of vehicle length; the main gear, on the other 
hand, may be placed either on the fuselage as a single unit, 
also as a fraction of fuselage length, or on the wing in multiple 
units. 

Wing 

The lifting planforms are assumed to be tapered, swept 
wings with straight leading and trailing edges. The pianform 
shape is trapezoidal as the root chord and tip chord are 
parallel. The wing is placed on the fuselage according to the 
location of the leading edge of its root chord, determined as a 
fraction of the fuselage length (Fig. 2). It is assumed that 
specified portions of the streamwise (aerodynamic) chord 
are required for controls and high lift devices, leaving the 
remainder for the structural wing box. The intersection of this 
structural box with the fuselage contours determines the 
location of the rectangular carrythrough structure. The width 
of the carrythrough structure is defined by the corresponding 
fuselage diameter. 
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For the transports in the present study, all the fuel is 
carried within the wing structure. An option is also available 
to carry the fuel entirely within Lhe fuselage, negating any 
bending relief in the wing. 

LOADS 

Fuselage 

Fuselage loading is determined on a station-by-station 
basis along the length of the vehicle. Three types of fuselage 
loads are considered — longitudinal acceleration, tank pressure, 
and bending moment. In the present study, all three load types 
are assumed to occur simultaneously to determine maximum 
compressive and tensile loads at the outer shell fibers at each 
station. 

Bending loads applied to the vehicle fuselage are obtained 
by simulating vehicle pitch-plane motion during a quasi-static 
pull-up maneuver; a landing; and movement over a runway 
bump. Simplified vehicle loading models are used where it is 
assumed that (1) fuselage lift forces (nominally zero for 
subsonic transports) are distributed uniformly over the 
fuselage plan area; (2) wing loading, determined 
independently, is transferred by a couple of vertical force and 
torque through the wing carry through structure; (3) fuselage 
weight is distributed uniformly over fuselage volume; 

(4) control surface forces and landing gear reactions are point 
loads; and (5) the propulsion system weight if mounted on 
the fuselage, is uniformly distributed. A factor of safety 
(nominally 1.5) is applied to each load case. The aircraft 
weight for each case is selected as a fraction of gross takeoff 
weight. All fuselage lift forces are assumed to be linear 
functions of angle of attack. Longitudinal bending moments 


are computed for each of the three loading cases and the 
envelope of the maximum values taken as the design loading 
condition. The bending moment computation is given in detail 
in Ref. 4 and will only be summarized here. 

Considering First the pull-up maneuver loading, the 
motion is assumed to be a quasi-static pitch-plane pull-up of 
given normal load factor n (nominally 2.5 for transport 
aircraft). The vehicle is trimmed with the appropriate control 
surface (a horizontal tail for all eight transport used for vali- 
dation in the present study), after which the angle of attack is 
calculated. 

Landing loads are developed as the aircraf t descends at a 
given vertical speed after which it impacts the ground; there- 
after the main and nose landing gears are assumed to exert a 
constant, or optionally a (1 - cos(cor)), force during ns stroke 
until the aircraft comes to rest. The vehicle weight is set equal 
to the nominal landing weight. Wing lift as a fraction of 
landing weight is specified, which reduces the effective load 
the landing gear carries. Likewise, the portion of total vehicle 
load the main gear carries is specified. No pitch-plane motion 
is coasidered during the landing. 

Runway bump loads are handled by inputting the bump 
load factor into the landing gear. Bump load factor is applied 
according to Ref. 8. This simulates the vehicle running over a 
bump during taxi. In a similar fashion to the landing, the wing 
lift as a fraction of gross takeoff weight is specified, as is the 
portion of effective load input through the main gear. No 
pitch-plane motion is considered during the bump. 
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Wing 

For the wing, only a quasi-static pull-up maneuver 
condition at load factor ft is considered for determining loads. 
At each span wise station along the quarter chord, from the 
wingtip to the wing-fuselage intersection, the lift load, center 
of pressure, inertia load, center of gravity, shear force, and 
bending moment are computed. For the inertia load, it is 
assumed that the fuel weight is distributed uniformly with 
respect to the wing volume. 

There is an option for either a trapezoidal or a Schrenk 
lift load distribution along the wingspan; the trapezoidal 
distribution represents a uniform lift over the wing area (which 
has a trapezoidal planform) while the Schrenk distribution is 
an average of the trapezoidal distribution with an elliptical 
distribution, where the lift is zero at the wingtip and maximum 
at the wing-fuselage intersection. Prandti has shown that a true 
elliptical lift load distribution will have a minimum induced 
drag, but a combination of the elliptical and trapezoidal 
distributions will give a better representation of actual aircraft 
loading. 9 


if x is behind it. Similarly, the total compressive stress 
resultant is 


n; = n xb + n Xa 


1 0, if not pressure stabilized 

N xpt if stabilized 


( 8 ) 


if x is ahead of the inlet entrance, and 

(0, if not pressure stabilized ] 




N 


X P' 


if stabilized 


(9) 


if x is behind it. These relations are based on the premise 
that acceleration loads never decrease stress resultants, but 
pressure loads may relieve stress, if pressure stabilization is 
chosen as an option. The stress resultant in the hoop direction 
is 


N y = rPgKp 


( 10 ) 


STRUCTURAL ANALYSIS 
Fuselage 

Weight estimating relationships are now developed for the 
load -carrying fuselage structure. In addition, the volume taken 
up by the fuselage structure is also determined. 

Considering first the circular shell, the stress resultants in 
the axial direction caused by longitudinal bending, axial 
acceleration, and pressure at a fuselage station x are 


Nx B /' 

(3) 

N = N * WS 

" X A p 

(4) 

AP. 

(5) 

* 

X 

“t) 

1! 


respectively, where r = D/2 is the fuselage radius, A = itr 2 is 
the fuselage cross-sectional area, and P = 2%r is the fuselage 
perimeter. In EQ (3), l' y = itr 3 is the moment of inertia of the 
shell divided by the shell thickness. In EQ (4), for the case of 
fuselage-mounted propulsion, W s is the portion of vehicle 
weight ahead of station x if x is ahead of the inlet entrance, or 
the portion of vehicle weight behind x if jc is behind the nozzle 
exit In EQ (5), P g is the limit gage pressure differential for the 
passenger compartment during cruise. The total tension stress 
resultant is then 


N + X =N XB+ N xp 


( 6 ) 


if x is ahead of the nozzle exit, and 

^x=N XB+ N Xp+ N XA 


(7) 


where K p accounts for the fact that not all of the shell material 
(for example, the core material in sandwich designs) is 
available for resisting hoop stress. 

The equivalent isolropic thicknesses of the shell are given 

by 


Wlu? 

II 

.2 

(ID 

II 

T 

(12) 

1 Sq - Kmg 1 mg 

(13) 


for designs limited by compressive yield strength (Fey), 
ultimate tensile strength (F lu ), and minimum gage, respec- 
tively. In EQ (13), t mg is a specified minimum material 
thickness and K mg is a parameter relating to l mg which 
depends on the shell geometry. 

A fourth thickness that must be considered is that for 
buckling critical designs, Fy „ , which will now be developed. 
The nominal vehicles of this study have integrally stiffened 
shells stabilized by ring frames. In the buckling analysis of 
these structures, the shell is analyzed as a wide column and the 
frames are sized by the Shanley criteria. 7 Expressions are 
derived for the equivalent isotropic thickness of the shell 
required to preclude buckling, i $ g , and for the smeared 
equivalent isotropic thickness of the ring frames required to 
preclude general instability. tp . The analysis will be restricted 
to the case of cylindrical shells. The major assumptions are 
that the structural shell behaves as an Euler beam and that all 
structural materials behave elastically. 

For the stiffened shell with frames concept, the common 
procedure of assuming the shell to be a wide column is 
adopted. If the frame spacing is defined as d and Young’s 
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modulus of the shell material is defined as E, the budding 
equation is then 


resultant N x on the cylinder. The buckling equation for these 
frameless sandwich shell concepts is 


n: 


dE 



(14) 




r 


Y" 


( 22 ) 


or, solving for 



Fuselage structural geometry concepts are presented in 
Table 1; values of the shell efficiency £ for the various 
structural concepts are given in Table 2. The structural shell 
geometries available are simply stiffened, Z-stiffened, and 
truss-core sandwich. We next size the frames to prevent 
general instability failure. The Shanley criterion is based on 
the premise that the frames act as elastic supports for the wide 
column; this criterion gives the smeared equivalent thickness 
of the frames as 



i K F id 3 Ep 


(16) 


where m is the buckling equation exponent. Or, solving for 
* S B 


( 





N 


-X 


rEz 


(23) 


This equation is based on small deflection theory, which 
seems reasonable for sandwich cylindrical shells, although it is 
known to be inaccurate for monocoque cylinders. Values of m 
and e may be found, for example in Refs. 10 and 1 1 for many 
shell geometries. Table 2 gives values for sandwich structural 
concepts available in PDCYL, numbers 8 and 9, both of which 
are truss-core sandwich. The quantities N x , r, and conse- 
quently t$ B , will vary with fuselage v station dimension x. 

At each fuselage station x, the shell must satisfy all failure 
criteria and meet all geometric constraints. Thus, the shell 
thickness is selected according to compression, tension, 
minimum gage, and buckling criteria, or 


where Cp is Shanley’s constant, Kp\ is a frame geometry 
parameter, and Ep is Young’s modulus for the frame material. 
(See Ref. 3 for a discussion of the applicability of this criterion 
and for a detailed derivation of the equations presented here.) 

If the structure is buckling critical, the total thickness is 

‘=‘s b +‘f b (17 > 

Minimizing t with respect to d results in 


*s = max ( f s c - l s T - h G • h B ) < 24 > 

If fj = i Sg , the structure is buckling critical and the 
equivalent isotropic thickness of the frames, tp, is computed 
from EQ (20). If i s > t Sg , the structure is not buckling critical 
at the optimum frame sizing and the frames are resized to 
make fc = h„. Specifically, a new frame spacing is computed 
from EQ (15) as 


4 1 

f ~ > 

nCp 

1 

8 

f 2 

2' P/H 

hi 

2 > 

ii 

a 

■ — 

l Kpx^EpE 3 J 

f 

) 

J 


(18) 


% 


‘F B 


= 7 1 (19) 

4 

= ( 20 ) 

4 


/ 



7 iCptE l 2 

KpiEp J 


( 21 ) 


where pp is the density of the frame material and p is the 
density of the shell material, so that the shell Is three times as 
heavy as the frames. 

Frameless sandwich shell concepts may also be used. For 
these concepts, it is assumed that the elliptical shell buckles 
at the load determined by the maximum compressive stress 


d = (25 

n; 

and this value is used in EQ (16) to determine ip 

The total thickness of the fuselage structure is then given 
by the summation of the smeared weights of the shell and the 
frames 


h = ' l S + ‘F < 26 ) 

The shell gage thickness may be computed from 
/" = t s / K mg . The ideal fuselage structural weight is obtained 
by summation over the vehicle length 

IV/ = 2!t^ (pi's, + PfQ; (27) 

where the quantities subscripted / depend on x. 

Since the preceding analysis gives only the ideal weight, 
W/ f the nonoptimum weight, W^q (including fasteners, 
cutouts, surface attachments, uniform gage penalties. 
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Table 1 Fuselage structural geometry concepts 


KCON sets 
concept number 

2 Simply stiffened shell, frames, sized for minimum weight in buckling 

3 Z-stiffened shell, frames, best buckling 

4 Z-stiffened shell, frames, buckling-minimum gage compromise 

5 Z-stiffened shell, frames, buckling-pressure compromise 

6 Truss-core sandwich, frames, best buckling 

8 Truss-core sandwich, no frames, best buckling 

9 Truss-core sandwich, no frames, buckling-minimum gage-pressure compromise 


Table 2 Fuselage structural geometry parameters 


Structural concept 
(KCON) 

m 

£ 

Kmg 

K p 

Kth 

2 

2 

0.656 

2.463 

2.463 

0.0 

3 

2 

0.911 

2.475 

2.475 

0.0 

4 

2 

0.760 

2.039 

1.835 

0.0 

5 

2 

0.760 

2.628 

1.576 

0.0 

6 

2 

0.605 

4.310 

3.965 

0.459 

8 

1.667 

0.4423 

4.820 

3.132 

0.405 

9 

1.667 

0.3615 

3.413 

3.413 

0.320 


manufacturing constraints, etc.) has yet to be determined. The 
method used will be explained in a later section. 


^ _ W BEND(y) 

?z s t 


(29) 


Wing 

Using the geometry and loads applied to the wing 
developed above, the structural dimensions and weight of the 
structural box may now be calculated. The wing structure is 
assumed to be a rectangular multi-web box beam with the 
webs running in the direction of the structural semispan. 
Reference 10 indicates that the critical instability mode for 
multi-web box beams is simultaneous buckling of the covers 
due to Local instability and of the webs due to flexure induced 
crushing. This reference gives the solidity (ratio of volume of 
structural material to total wing box volume) of the least 
weight multi-web box beams as 


1 = 



) 


(28) 


where W bend is the weight of bending material per unit span 
and p is the material density. W bend is computed from 
EQS (28) and (29). The weight per unit span of the shear 
material is 


W SHEAR 



(30) 


where is the applied shear load and a$ is the allowable 
shear stress. The optimum web spacing is computed from 2 


dw =' 


(»- 2 e c ) 

M ] 

(l-ec)V 2E tV 

Us' 2£ J 


2e c -l 

2e c 


3 

t c Ze c 


2 


where e and e depend on the cover and web geometries (3 1 ) 

(Table 3), M is the applied moment, t is the thickness, E is the 
elastic modulus, and Z 5 is obtained from Ref. 10. The solidity 
is therefore 
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where subscripts W and C refer to webs and covers, 
respectively. The equivalent isotropic thicknesses of the 
covers and webs are 


Table 3 Wing structural coefficients and exponents 


Unstiffened 

Unstiffened 

Unstiffened 

Truss 
Trass 
Trass 


rc "M 


Webs 

E 

e 

E 

e c 

Ew 


Truss 

Unflanged 

Z-stiffened 

2.25 

2.21 

2.05 

0.556 

0.556 

0.556 

3.62 

3.62 

3.62 

3 

3 

3 

0.605 

0.656 

0.911 

1.000 

1.000 

1.000 

Truss 

Unflanged 

Z-stiffened 

2.44 

2.40 

2.25 

0.600 

0.600 

0.600 

1.108 

1.108 

1.108 

2 

2 

2 

0.605 

0.656 

0.911 

0.546 

0.546 

0.546 


(32) 




M 


K e c) 


Z s i 2 E 


I' ecdw Vr 


t ) 


2 

lew 


— (33) 


respectively, and the gage thicknesses are 


~ ^gc* c 

(34) 

= K gw i W 

(35) 


Values of e, e, to E C. ew, ““f/fc •* found '" 
Table 3 for various structural concepts. If the wing s 
semispan is divided into N equallength segments, the total 
ideal weight of the wing box structure is 


W B OX = ^f-^( W BEND, + W SHEAR , ) 
i=l 


(36) 


The wing carrythrough structure consists of torsion 
material in addition to bending and shear material. The torsion 
material is required to resist the twist induced due to the sweep 
of the wing. The bending material is computed in a similar 
manner as that of the box except that only the longitudinal 
component of the bending moment contributes. Letting 
f0 = r(y = 0) and Mq - M(y = 0), 




( M 0 cof(hs) 


The weight of the bending material is then 
w B end c = P*C C SR l 0»C 


(37) 


(38) 


where w c is the width of the carrythrough structure. (When 
the wing-fuselage intersection occurs entirely within the 
cylindrical midsection. as is the case with all eight transport 


0.407 

0.505 

0.405 


used for validation in the present study, wc = D.) The quanti- 
ties dw- rw. and fc are computed in the same manner as for die 
box. The weight of the shear material is 


FSo 

WSHEAR C = P ^7 WC 

where Fj n = F$(0). u 

The torque on the carrythrough structure is 

7=M 0 sin(As) 

and the weight of the torsion material is then 

_ pr(f 0 + c 5 ^)H' C 

WTORSIONc - (q c sr o s 


( 39 ) 


( 40 ) 


( 41 ) 


Finally, the ideal weight of the carrythrough structure is 
computed from a summation of the bending shear and torsion 
material, or 

W C = W b END c + W SHEARc + ORS!ON c (42) 

As in the case of the fuselage structural weight, non- 
optimum weight must be added to the ideal weight to obtain 
the true wing structural weight. The method used will be 
discussed below. 

REGRESSION ANALYSIS 

Overview 

Using fuselage and wing weight statements of eight 
subsonic transports, a relation between 
bearing structure weights obtained through PDCYL and the 
actual load-bearing structure weights, primary structure 
weights, and total weights is determined using statistical 
analysis techniques. A basic application which is first 
described is linear regression, wherein the esuma ed wights 
of the aircraft are related to the weights calculated by PDCYL 
with a straight line, y = mi + b, where y is the value of the 
estimated weight, m is the slope of the line, x is the v u 
obtained through PDCYL. and b is the y-intercept. This line is 
termed a regression line, and is found by using the method of 
least squares, in which the sum of the squares ot the residual 
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errors between actual data points and the corresponding points 
on t be regression line is minimized. Effectively, a straight line 
is drawn through a set of ordered pairs of data (in this case 
eight weights obtained through PDCYL and the corresponding 
actual weights) so that the aggregate deviation of the actual 
weights above or below this line is minimized. The estimated 
weight is therefore dependent upon the independent PDCYL 
weight 

Of key importance is the degree of accuracy to which the 
prediction techniques are able to estimate actual aircraft 
weight A measure of this accuracy, the correlation coefficient 
denoted /?, represents the reduction in residual error due to the 
regression technique. R is defined as 


/? = 


E t -E r 

E r 


(43) 


PDCYL analysis only predicts the load-carrying structure of 
the aircraft components, a correlation between the predicted 
weight and the actual load-carrying structural weight and 
primary weight as well as the total weight of the fuselage, was 
made. 

Structural weight consists of all load-carrying members 
including bulkheads and frames, minor frames, covering, 
covering stiffeners, and longerons. For the linear curve-fit the 
resulting regression equation is 

^actual = L3503 W calc R = 0.9946 (46) 

This shows that the nonoptimum factor for fuselage 
structure is 1.3503; in other words, the calculated weight must 
be increased by about 35 percent to get the actual structural 
weight For the alternative power- intercept curve fitting 
analysis, the resulting load-carrying regression equation is 


where E t and E r refer to the residual errors associated with the 
regression before and after analysis is performed, respectively. 
A value of R = 1 denotes a perfect fit of the data with the 
regression line. Conversely, a value of R = 0 denotes no 
improvement in the data fit due to regression analysis. 

There are two basic forms of equations which are 
implemented in this study. The first is of the form 

y e « = "Ucalc ( 44 ) 

The second general form is 

yest = ^calc (45) 

Fuselage 


w 'actual = 1.1304W‘2 179 R = 0.9946 (47) 

To use either of these equations. to estimate total fuselage 
weight, nonstructural weight items must be estimated inde- 
pendently and added to the structural weight. 

Primary weight consists of all load-carrying members as 
well as any secondary structural items such as joints fasteners, 
keel beam, fail-safe straps, flooring, flooring structural 
supplies, and pressure web. It also includes the lavatory 
structure, galley support, partitions, shear ties, tie rods, 
structural firewall, torque boxes, and attachment fittings. The 
linear curve fit for this weight yields the following primary 
regression equation 

'actual =1.8872^ R = 0.9917 (48) 


The analysis above is used to develop a relationship 
between weight calculated by PDCYL and actual wing and 
fuselage weights. The data were obtained from detailed weight 
breakdowns of eight transport aircraft 12 * 16 and are shown in 
Table 4 for the fuselage. Because the theory used in the 


The primary power- intercept regression equation is 

^actual = 1.6399 w!^ 141 R = 0.9917 (49 ) 


Table 4 Fuselage weight breakdowns for eight transport aircraft 


Weight, lb 

Aircraft 

PDCYL 

Load-carrying structure 

Primary structure 

Total structure 

B-720 

6545 

9013 

13336 

19383 

B-727 

5888 

8790 

12424 

17586 

B-737 

3428 

5089 

7435 

11831 

B-747 

28039 

39936 

55207 

72659 

DC-8 

9527 

13312 

18584 

24886 

MD-11 

20915 

25970 

34999 

54936 

MD-83 

7443 

9410 

11880 

16432 

L-1011 

21608 

28352 

41804 

52329 
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Tbe total fuselage weight accounts for all members of the 
body, i iKt'otiTig the structural weight and primary weight. It 
does not i iyiwde passenger accommodations, such as seats, 
lavatories, kitchens, stowage, and lighting; the electrical 
system; flig ht and navigation systems; alighting gear, fuel and 
propulsion systems; hydraulic and pneumatic systems; die 
communication system; cargo accommodations; flight deck 
accommodations; air conditioning equipment; the auxiliary 
power system; and emergency systems. Linear regression 
results in the following total fuselage weight equation 

^actual = 2.5686W caic R = 0.9944 (50) 

This shows that tbe nonoptimum factor for the total 
fuselage weight is 2.5686; in other words, the fuselage 
structure weight estimated by PDCYL must be increased by 
about 157 percent to get the actual total fuselage weight This 
nonoptimum factor is used to compare fuselage structure 
weight estimates from PDCYL with total fuselage weight 
estimates from tbe Sanders and the Air Force equations used 
by ACSYNT. 

The total fuselage weight power-intercept regression 
equation is 


W^=3.9089W™ 578 


R = 0.9949 


(51) 


Plots of actual fuselage component weight versus 
PDCYL -calculated weight, as well as the corresponding linear 
regressions, are shown in Figs. 3-5. 


Wing 

Tbe same analysis was performed on the wing weight for 
the sample aircraft and is shown in Table 5. The wing box, 
or load-carrying structure, consists of spar caps, interspar 
coverings, spanwise stiffeners, spar webs, spar stiffeners, and 
interspar ribs. The wing box linear regression equation is 

W, actual = 0.9843 W cak R = 0.9898 (52) 

so that the nonoptimum factor is 0.9843. Power-intercept 
regression results in 

W actua i=\3342W^ 01 R = 0.9902 (53) 

Wing primary structural weight includes all wing box 
items in addition to auxiliary spar caps and spar webs, joints 
and fasteners, landing gear support beam, leading and trailing 
edges, tips, structural firewall bulkheads, jacket fittings, 
terminal fittings, and attachments. Linear regression results in 

W actual = 1-3442 W caic * R = 0.9958 (54) 

Power-intercept regression yields 

^ actual = 2.1926W™ 534 R = 0.9969 (55) 


o 
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DC-8 


MD-11 

MD-83 


Table 


5. Wing weight breakdowns for eight transport aircraft 



22080 


33617 

6953 


35157 

8720 

28355 


27924 


47614 

11553 

36101 


35330 


62985 

15839 

46233 


' The total wing weight includes wing box and primary 
weight items in addition to high-lift devices, control surfaces, 
and access items. It does not include the propulsion system, 
fuel system, and thrust reversers; the electrical system, 
alighting gear; hydraulic and pneumatic systems; ami-tang 
Sesfand emergency systems. The resulung total weight 

linear regression equation is 

^actual ~ l-7372W fa / r « = 0.9925 (56) 


This shows that the nonoptimum factor for the to wing 
weight is 1.7372; in other words, the wing box weight esti- 
nated by PDCYL must be increased by about 74 percent 
get the actual total wing weight This nonopumum factoMS 
used to compare wing box weight estimates from PDCYL with 
total wing weight estimates from the Sanders and the Air 

Force equations used by ACSYNT. 

The power-intercept equation for total wing weig is 

3.7464W£T* 09946 (5?) 


Plots of actual wing component weight versus PDCYL - 
calculated weight as well as the corresponding linear 
regressions, are shown in Figs. 6-8. 


Discussion 

Both fuselage and wing weight linear and power regres- 
sions give excellent correlation with the respective weights 
of existing aircraft, as evidenced by the high values ot the 
correlation coefficient, R. It should be noted that even thoug 
Oie power-based regressions give correlations equal to or 
better than the linear regressions their factors may va D 
distinctly from the linear cases. This is due to their po 

n<H ^ecau^ estimates of non- load-bearing primary structure 
are generally not available at the conceptual design stage and 
because nonprimary structure is probably not well esumated 
nonoptimum factor, EQS (48) and (54) are recommended 
for estimadng the primary structural weights of the respecuve 
rran snort fuselage and wing structures (Figs. 4 and /). 

may be made be, ween we,gh, esumams 

from weigh, esumaUng relauo.ish.ps currency used by 
ACSYNT PDCYL output, and actual aircraft componen 
weights. Figure 9(a) shows a comparison between fttsdage 
weight estimated from the Sanders equauon. the Air Force 
equation, and PDCYL output with the actual tuselage weight 
of the 747-21P. Figure 9(b) shows a similar comparison 
wing weight. 
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Actual Primary Wing Weight Actual Wing Box Weight 
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Fig. 6 Wing load-carrying structure and linear regression. 
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Fig. 7 Wing primary structure and linear regression. 
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Fig. 9(b) Wing weight estimation comparison for 747-21 P. 
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ADDING COMPOSITES TO PDCYL 

Light weight materials such as fiber reinforced plastics (composites) and 
bonded honeycomb sandwiches have become more and more common in 
airplanes in the last two decades (1). Designers value the unique 
properties of these materials, particularly their high stiffness to weight 
ratios. They must, however, balance these assets against the additional 
cost of these materials and their manufacture. To aid designers with this 
analysis, a composites subroutine has been added to the PDCYL structures 
weight calculation code. This subroutine sizes the thickness of a 
particular composite necessary to withstand the required aircraft loads, 
and provides this information to PDCYL which calculates the resultant 
weight of the aircraft. 


TYPICAL AIRCRAFT COMPOSITES 

The selection and use of composites on transport aircraft is an evolving 
process. A variety of composites have been tested in both military and 
commercial aircraft in the last 25 years (1). These composites typically 
consist of a strong, stiff fiber such as glass, graphite or kevlar, and a 
protective, adhering, inexpensive plastic matrix such as polyester or 
epoxy. 

Glass fibers embedded in a polyester matrix have been the dominate 
composite for military and civil aircraft in the past. Currently, the 
aircraft industry prefers the stiffer and higher temperature composites 
made from carbon fiber in an epoxy matrix. However, the grade of carbon 
fiber and epoxy seems to change from year to year and from airplane 
manufacturer to airplane manufacturer. The current favored carbon fibers 
are AS4 (Hercules/Hexcel), IM6 and IM7 (Hercules/Hexcel) . The AS4 is an 
economical, high-strength carbon fiber and the IM6 &7 are high-modulus 
expensive fibers. These three carbon fibers have been used on military 



aircraft and in research, but are not on commercial vehicles. The T-800 
fiber (a Toray equivalent to IM7) has recently been used in some 
commercial applications (1-6). 

Epoxies, particularly the 350°F curing systems, are the least expensive 
high temperature options for matrix materials. Several epoxy systems 
have been developed and tested for use with specific fibers. There is a 
current trend to use rubber modified epoxies such as 8552 and 3900 to 
increase the toughness of the composite system and its resistance to 
impact. Fiber-resin combinations currently in use by airplane designers 
and researchers are: 

AS4/938 (ICI Fiberite) -Boeing Advanced Composites Program Door 
Panel(2) 

AS4/8552 (Hexcel/Hercules), -Boeing Adv. Comp Fuselage (6-7) 
AS4/8551 (Hexcel/Hercules) (6) 

AS4/3501-6 (Hercules) -McDonnell Douglas Adv. Technology 
Composite Wing program (8) 

AS4/3502 (Hercules) Military Aircraft (6) 

COMPOSITE STRUCTURAL ANALYSIS 

Composite materials were originally added to the options in the 
PDCYL program in 1995. This was done by simulating these materials by 
homogeneous structures with uniform mechanical properties (strength 
and modulus of elasticity) in every direction. This approach limits the 
code to evaluation of only the simplest and weakest type of composites 
called random mati. Random mat composites are made by stacking the 
reinforcing fiber in all direction throughout the thickness of the 
material. In this type of composite the elastic properties and strength of 
the layup are roughly the same in every direction but the fiber density and 
reinforcement is low in any specific direction. 

Random mat composites are not favored by aircraft designers because of 
their low strength to weight ratios. The preferred type of composite for 
these applications are ones in which the properties of the material are 
customized to meet the specific directions and magnitudes of the 
structural loads. This yields the minimum weight composite for the job. 

To accomplish this, composite designers specify a layup pattern for a 
composite laminate relative to a major axis of loading. 


i See Appendix A for definition of composite terms 




A typical composite laminate is made of a stack of 4-16 plys. A ply is a 
single layer of parallel reinforcing fibers embedded in a partially cured 
matrix of plastic. The location of each ply in the stack is defined 
relative to the angle its fibers makes with a major axis, such as the x- 
axis. For instance, a 0/90/90/45/0 layup is one in which the fibers of the 

outer and inner layers are parallel to the x-axis, the next two plies have 

fibers perpendicular to this axis and the fibers of the third layer are at an 

angle 45° clockwise to the x axis. This type of composite would have 

reinforcing fibers to sustain tensile and compressive loads in the x and y 
directions but would be weakest in the 45° direction. Composites walls 
for structural parts such as aircraft are often made from stacks of these 
laminates. 

Analysis of a multilayer stack is more complex than that of homogeneous 
materials such as aluminum or random mat and requires the use of a 
macromechanics approach to determine elastic properties and strength. 
The macromechanics approach used in the COMPOS part of the PDCYL code 
is that presented in most textbooks on composite design (9-11) . In this 
approach the stiffness of a particular laminate is calculated by summing 
the contributions of each layer (ply) in the stack to the stiffness of the 
laminate in a particular direction. The composite stiffness in each major 
direction is then used to calculate the net strain of the composite in that 
direction due to the applied loads. From the net strain, the strain on each 
layer (ply) parallel and transverse to its fiber is derived. The resulting 
strains are then compared to the failure strains of the ply material and 
from this the potential for the failure of the stack is determined. The 
details of implementing this approach in PDCYL are described in the 
following section describing the COMPOS (composites) code addition. 

COMPOS CODE 

COMPOS is a section of code which has been added to PDCYL program to 
calculate the minimum laminate thickness required to withstand the 
forces imposed at each section of the airplane. 

Assumptions within COMPOS 

•The laminate is symmetric and orthotropic. (This type of layup is 
commonly used in aircraft design to minimize warpage of the layup). 

• Every ply in the stack is composed of the same resin- fiber material. 

• The stack is a minimum of 3 plys. (A ply is usually .003-. 007 inches 
thick depending on the material.) 



• The modulus of the material is the same in compression and tension, (if 
the compression modulus is different than its tensile modulus, the 
smaller of the two values is selected for all calculations.) 

• Failure of the composite laminate occurs when any single ply fails. 

• Failure of a ply occurs when it reaches the maximum strain transverse 
or parallel to the fiber direction in tension, compression or shear (11) 
Maximum strain theory is invoked in this analysis because it is currently 
believed to be the most predictive failure theory for composites (3,4,8) . 
•The minimum gage thickness for the composite material is assumed to 
be the thickness of the initial laminate (a stack of plies). 

•All loads are applied in the plane of the ply. This means that there are 
no z direction loads in tension, compression or shear. 

• The buckling equations used in PDCYL to analyze the frames and 
stringers made from homogeneous materials apply to these heterogeneous 
materials. For buckling analysis the modulus of the laminate in the 
direction of load is used. This is a very course assumption and maybe 
somewhat optimistic for quasi isotropic composites manufactured with 
adhesive joints but seems highly unlikely for symmetric orthotropic 
laminates with heterogeneous properties. However, buckling analysis of 
complex composites structures is still in the developmental stage. (12) 

Calculation Procedure 

•Calculations for Compressive and Tensile Loads 
Once the maximum tensile and/or compressive loads per unit width (Nx 
and Ny) at any given aircraft section are determined in the PDCYL code, 
they are transferred to the COMPOS subroutine. The effect of these 
normalized forces on the composite laminate strain is calculated using 
the following relationship for an orthotropic symmetric laminate (9) : 


[N] = [A] x [e°] 


( 1 ) 


Where: 

[N] = Matrix of forces on the composite section (Nx, Ny and Nxy) 

[A] = Stiffness matrix of the composite 

[e°] = strain matrix of the composite (ex, ey, exy) 


The components of the stiffness matrix [A] are determined in the code 
through the following relationship (9): 


Ajj=X n K =i (QBy)k(hk) 


(2) 



Where: 

QBjj = component of each ply’s stiffness in the i and j’s directions 

h - thickness of k ply 

k= ply number in the laminate 

The stiffness contributions, QB values, of each ply are determined from 
the initial ply properties, El, E2, v 12 and the ply angles, 0, specified by 
the user in the input file for a particular laminate construction. (Here, the 
“1” direction is taken parallel to the fiber and the "2” direction 
transverse to the fiber). 

Once the average laminate strain is determined from equation (1), this 
strain is then transferred to each ply and transformed into strain parallel 
and transverse to each fiber as well as shear strain. These strains are 
then divided by the mating failure strains for the material (supplied by 
the user in the input file) to determine the R value of the layup. 

Rjj= alley /e jj 

Where: 

alley-allowable components of strain in principle ply direction 
e jj = components of strain in principle ply directions 

If the R value for all plys in all the principle directions is more than 1, 
the laminate thickness is adequate to support the load and is left 
unaltered. If R is less than one on any ply in any of the principle 
directions, the thickness of the laminate is increased by giving it the 
value of it initial thickness divided by R. 

Calculations for Buckling 

PDCYL currently determines critical buckling loads from the modulus of 
elasticity of the material. COMPOS calculates the modulus of the 
laminate in the direction parallel to the buckling force and passes this 
value back to PDCYL. As mentioned in the assumptions portion of this 
report, the buckling calculation of PDCYL may not be valid for composites 
as they were developed for isotropic materials. Little research has been 
done on composites in buckling so the authors advise caution in 
interpreting this result particularly with non-isotropic layups. 


NON-OPTIMUM FACTORS 



Unfortunately, few all composite planes have been built so it is 
difficult to find planes to use as checks for the composite section of the 
code (8). The all composite planes listed in the literature (8) are: 

• Windecker Eagle in 1967 which was glass fiber reinforced 

• Learfan in 1981 which used glass, carbon and kevlar fibers 

• Piaggio Avanti in 1986 with carbon fiber parts 

• Beech Starship in 1986 with carbon fiber 

• Grob GF-200- all composite 

• Slingsby T-3A Firefly -all composite 

A literature search and personal interviews failed to turn up much 
information directly useful in determining non-optimum factors. (These 
factors are used to multiply the results of theoretical calculations to get 
weights of practical structures.) 

One reference was found which had this type of data (12). In this 
reference, a theoretical analysis gave 8640 pounds as the weight of a 
composite wing box whereas the actual wing was estimated to weigh 
11,284 pounds giving a non-optimum factor of 1.306. Using the non- 
optimum factors for aluminum structures (13) this number can be used to 
estimate non optimum factors for carbon fiber-epoxy structures. If It 
is assumed that the non-optimum factors for the fuselage primary 
structure increase in the same proportion as wing structure relative to 
aluminum, and that the increments for secondary structure and non- 
structural are the same for graphite-epoxy composites and aluminum, 
then the following non-optimum factors for the composite result: 



Primary 

Primary ASecondary 

Total 


Structure 

Structure 

Assembly 

Fuselage 

1.792 

2.329 

3.010 

Wing 

1.306 

1.666 

2.059 


There are many composite components in commercial and military 
structure as well as some from research on advanced composites. It may 
be possible to compare these components to predictions of the code. 
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APPENDIX 1 

COMPOSITES TERMINOLOGY 


Random Mat - equal fibers in every direction 

Balanced - equal fibers in orthotropic directions yield a composite with 
identical properties in 2 principal directions. 

Symmetric -A symmetric laminate is one in which for each ply above the 
center of the stack there is and identical one at an equal distance below 
the center. For instance, a 0/-45/90/90/-45/0 is a symmetric layup 
but a 0/-45/90/ 0/-45/90 is not. 

Quasi-lsotropic - Layups which are designed to have only two independent 
elastic constants, the modulus of elasticity and Poisson’s ratio. These 
materials have the same values in every inplane direction. To meet this 
criteria fiber (ply) layups must have the following conditions: 

•Total number of plies must be 3 or more 

•Individual plies must have identical stiffness [Q] matrices and 
thickness 

•Layers must be oriented at “equal” angles (if total number of 
layers is n, than each layer is pi/n relative to the next). If the 
laminate is constructed from several groups of laminates,, the 
condition must be satisfied for each laminate group 
Typical laminates which satisfy these rules : [0/60/-60], [0/45/- 
45/90] 



Appendix C 


Weights of ASA 2150 With Aluminum and Graphite/Epoxy Fuselage 
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STRUCTURAL CKF FSK EFF CK PG CF CTHIC SAFEFAC DEFL 

PARAMETERS 5.2400 0.04504 0.76000 2.0390 11.250 0.6250E-04 0.000 1.500 0.000 
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Approximate Methods of Trajectory Optimization 
Mark D. Ardema 


Introduction 

Application of optimal control theory in the form of the maximum principle to aircraft 

trajectory optimization problems generally results in a two-point (2PBVP). The order of 

A 

this problem is double the number of state variables and the equations are always “half 
unstable.” Many schemes have been developed to numerically solve this difficult class of 
problem, but all are unsuitable in a vehicle synthesis code. Not only are they 
computationally expensive, but they are non-robust and not user-friendly. 

What is needed in a vehicle synthesis code is a method that optimizes the trajectory in 
one pass, that is as an integral part of the trajectory integration. The method must also be 
robust and it should be easy to use and to interpret physically. The key to achieving this is 
to use judicious approximations to reduce the functional optimization problem to a 
function one. 

In this report, two approximation techniques are reviewed and developed. The first is 
the use of the Energy State Approximation (ESA). This well-known technique substitutes 
the total mechanical energy for the speed as a state variable, and then neglects the altitude 
and flight path dynamics relative to the energy dynamics. The second technique is the use 
of Singular Perturbation Theory (SPT) to time-scale decouple equations of motion. These 
two techniques are related, and in fact, the ESA may be viewed as an example of SPT 
methods. 



Trajectory O ptimization- the Maximum Principle 


The equations of motion of aircraft flight, no matter what the assumptions (see 
Appendix A), are of state equation form: 

x = L{2LU) 


where Xe^ a is the state and U sUCSR" is the control vector. Suitable boundary 
conditions on the state vector components are prescribed. It is desired to find the 
components of U along the trajectory so that 

J = £' fo(X, U)dt 

is minimized. It is assumed that time is free. The necessary conditions for optimal control 
are provided by the Maximum Principle (MP). 


Theorem (the Maximum Principle): Introduce the variational Hamiltonian function 

H = -i + jVi 

/=l 

where the components of the adjoint vector, X t - , satisfy the differential equations 



i = 1, ., n 


Then, if U_ is an optimal control. 


(a) U_ = arg max H 
U_ eU 


(b) H = 0 

(c) Transversality conditions (“natural” boundary conditions on the A t ) hold. 
Thus, we must solve a 2jj dimension 2PBVP in the states and adjoints; exactly n 


2 


boundary conditions are provided at / = O and the other half at / = tf (due to the 
travsversality conditions). The equations are unstable in the sense that if they are 
linearized about a nominal trajectory, one-half of the eigenvalues will have positive real 
parts and the other half negative (unless they are zero). 

A pproxim ation Te chniques 

Methods of reducing the 2PBVP to a simpler problem will now be developed. These 
methods focus on order reduction and are motivated by two simple observations. 

First, we note that if all componets of / , except possibly 4 and the function 4 are 
independent of a specific state variable, say X* and the final value of X, is not specified, 
then the corresponding adjoint is always identically zero and the state equation X , = f i 

drops out of the problem (decouples from the other states). To see this, consider the i* 
adjoint equation and its travsversality condition: 

-i ,—•££*, M'rh o 

d X i 

The only solution to this problem for any finite value of d f t jdX t , is X i s 0. 

Second, we note that if there is only one state equation, then the necessary conditions 
can be used to eliminate the adjoint variable and thus the problem reduces from a 
functional optimization problem to a function one. To see this, consider 


X = f(X,U) 



X = ?JL-il±- 

d X d X 

Applying the MP (assuming for the moment unbounded optimal control exists) 


H = -fo + X f = 0 


d H d fo , d f . 

d U d U d U 

Eliminating X -gives 

d fo r dff A 

/ + fo - 0 

d U d U 

for the optimal control. Alternatively, a direct approach may be used: 

■'■ = /£<* = 

Thus {fo/f ) is to be minimized with respect to U at constant X ; this leads directly to the 
equation for optimal control derived just above from the MP. 

SPT provides an organized, mathematical way to view order reduction of differential 
equations. Consider the initial value system 


X = f(X, Y ) X(0) = X, 

sr=g(x, r) (0)=n 

where e is a “small” parameter. 

Since e is small, an approximate system may be expected to be 
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x,=f(x„r,) 


0 = g(X„Y r ) 


It can be proved that under certain conditions, the solution of this problem is a good 
approximation to the solution of the original problem, except near t = 0 because the 
boundary condition y(0) = Y 0 will be generally violated these. 

The problem is that ^undergoes a rapid transition from its boundary condition to the 
approximate solution at t = 0 . To analyze this motion, the time scale is stretched by 
T -tie. The resulting equations are called the boundary layer equations 


dX 

dT 


= ef{X,Y) 


dY 

dT 


S(X,Y) 


Setting e = 0 to approximate these equations 


= 0 => X = const = X 0 

dT 0 


-jjr = g(X t ,Y t ) 

The solution to this equation approximates the desired solution near t = 0 . There are 
matching techniques to combine these two solutions to give an over-all approximation, if 
desired. The key observation is that a second order system has been replaced by two first 
order systems, and each of them reduces to a function optimization problem. 
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The SPT provides a convenient way to look at the energy state approximation. Define 
the aircraft energy per unit weight by 

E = h + — V 2 

2g 

Differentiate and use the state equations in the Appendix 


g Mg 


Where T v is the component of thrust along _V and P is the specific excess power. Note 
that this equation is valid for all three sets of equations given in the Appendix. 

Now replace V by E as state variable and use the observation that h and y are 
capable of rapid change relative to E . This motivates writing 

E = ? 


s h —• 


£/=— 

Setting e - 0 then gives an order reduction of the equations of motion by two. We will 
use this approximation, the ESA, throughout. 

This approximation has a long history of successful application in a wide variety of flight 
trajectory problems. The main drawback is that the variables h and y may now jump 
instantaneously at points along the trajectory, as well as at the boundaries. These jumps 
could be accounted for by boundary layer analysis, but this is not done in this report. 
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ESA equations of motion for the three cases of interest are given in the Appendix, 
these equations will now be used as the basis for discussing specific trajectory 
optimization problems. 

Minimum Time/Fuel to Climb 

Starting from equations (1)' 


m = -CT zn(0) = m 0 

X - V 

E = P E(0) = E 0 ,E(t f ) = E f 

L = mg 

/ = \{K, +K 2 CT)dt 

Here, the system functions and boundary conditions do not depend on X and thus the 
equation X - V drops out of the problem: 

m = -CT m( 0) = m 0 

E = P E(0) = E 0 ,E(t f ) = E f 


/ = J (AT, +K 2 CT)dt 

with P evaluated at L- mg . The system functions and boundary conditions now 
depend on both E and M and hence neither state equation uncouples. Thus the MP must 
be applied and a 2PBVP solved. 
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To reduce the problem to one of function optimization, SPT is used to further reduce 
the system. 

m - -CT 


eE = P 

Setting e = 0 gives a single state equation 
m = -CT /n(0) = 

T - D, L-W 

The optimization problem is now (h and E are controls) 

J = (K, +K 2 CT) t f 

With the obvious trivial solution t f = 0 . Also, because the system functions do not 
depend on m and m(^ / ) is free, A M = 0 . 

The boundary layer system for this problem with s = 0 is simply 

E = P 

with m = const , so that 

J = j(Kx + K,CT)y. 
and the solution reduces to 


maxf p N 
h [k, +K 2 CT; 


E = const 
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assuming that 



is positive and 2: is monotonic. This is the well-known 


energy climb path. 


From now on we will assume “slowly varying” mass, that is that m is on a slower time 
scale than E and thus its state equation may be ignored. It is also assumed that the throttle 
is fixed. 


Fixed Ranee 


This problem is the same except that range is fixed 


X = V 


x(0) = x„x(t r ) = x r 


E = P 


E(0) = E„E(t r ) = E r 


S = f(jC, +K,CT)dt 

Here, X x = const ^ 0 so that the X = V state equation does not uncouple, and we have 
a 2PBVP. To effect system order reduction, SPT is used. 

X = V 

e E = P 

H - -K x - K 2 CT + X x V + X E P 

• d H . ■ , 

L= = 0 => X y = const 

x d X 
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sX B — 


a h a(cf) a v ar 

d E 2 d E x d E ' d E 


Setting £ =0a problem with a single state is obtained 


X - V 


T=D 


H = -AT, -K 2 CT + X x V 


L=W 


Applying the MP: 

Min (K, + K 2 CT\ _ K, +K 2 C c T c 

h,E V V ) l-w = V c 


K, +K,C c T c 

x ~ V 
Y c 

This defines a cruise point, characterized by C c , T c , and in the flight envelope. By 
proper selection of Ki and K 2 , this point can be made to closely approximate minimum 
direct-operating-cost cruise. 

For minimum time (AT, = 1, K 2 = 0) , the optimum cruise point is given by 


Max 

h,E 

For minimum fuel consumption (AT, =0, K 2 = l) it is 



Max 

h,E 


( V'] Max 


Icf) ~ h,E 

1 c J 


which is the classic Brequet cruise point. 
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The boundary layer with s = 0 is 


E = P 


H = -K x -X 2 CT + X x V + X B P 
so that the optimal climb flight path is given by 


Maxf P > 

h [Kt Tk^CT^TxV; 


E- cotut 


with X x as given above. 


Maximum Turning With No Thrust 


We start with (2)' with T=0. 


X = V cos x 


Y = V sin x 


E = ? = - 


VD 

mg 


L sin a 

X = 

mV 


L cos <f> = mg 


n 



with 


„ f . . r L sin 6 


In this case the X, Kand equations all uncouple. Changing to more convenient 
variables: 


E = -v(b + Qo 2 ) 

J - 0) fdt 

where 


<S) = 


tan# 

tan# v 


, - 1 <, co £ +1 


<f> M = sec" l [min(<7 t limit, load factor limit)] 


Thus 


f _ g tan ^ 
V 




<of 


v(B + Ca> 2 ) 


dE 


and the optimal controls are given by 


Max 

b,ca 


f 


cof 


v(b+cw 2 ); 


Iff* coast 
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which leads to 



=> a) = 




(assuming right hand turn) 


d_ 

dh 


() = 0 


(v— + £) (£ + Ca 2 )-fv(— + — <v 2 )= 0 
V. d h VJ ' > \d h dh ) 


where 

B Dq±Dlo c Dj^_ 
Mg ’ Mg 3 


so that 




The search for the optimum h probably should be done numerically. 

Next, consider the same problem but using (3)' . Now, the X equation uncouples 
but the Y and x equations do not. The coupling is of two types. First, through the 
Coriolis terms, which are relatively small and can be ignored. Second, through the 
centripetal terms, which are large at the start of descent trajectories from orbit. 
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There are two ways to deal with this problem. First, the Coriolis and centripetal 
terms are ignored. This is justified because what is really sought is turning ability due to 
banking and these terms mask this. Second, the E and x terms may be decoupled using 
SPT. 

Maximum Cross-Ranee 

Next consider, using (2)' 


f = -j Ydt =-jVsin X dt 


X - V cos^ 

Y=V sin* 

E = P 

L sin 6 

X = 77~ 

wV 

Lcosfi = mg 

As before, the Xand Y equations uncouple but now the x equation does not. To 
reduce this to a function optimization problem, further time-scale separation is required. 
Putting x °n a slower time scale than E gives the solution V = 0 and X x = 0. 

t 

Using equations (3) results in the same problems as for maximum turning. 
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Appendix 

Equations of Motion 


The following are the aircraft point-mass equations of motion under various 
approximations. 

(1) Flight in a vertical plane over a flat, non-rotating earth; no winds aloft and 
thrust aligned with velocity. 

m - -CT 

X = V cos y 

h - Vsiny 

y - H ~ mgsxny 

m 

. _ L- mg cosy 
7 ~ ~mV 

(2) 3-D flight, otherwise the same as (1). 

m - -CT 
X = Vcosy sin/ 

Y = V cos^ cos/ 
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h = sin?' 


V = 


T-_D-mgmiy 

1 37 


L sin^ 
mV cosy 


. Lcosi- mg cosy 

y = 

1 mV 

(3) 3-D flight over a spherical, rotating earth; no winds aloft, thrust not aligned 
with velocity, terms in the square of the earth rotation ignored. 

T 

m = 

Ss^sp 


V = 


T cos/3 cos (a+£)-D 

m 


-gsiny 


(T cos/3 sin(a + £) + L 
l mV 


e V 

cost-2-cosy + — cosy + 2n> cos cos 


f T cos/3 sin(a + £)+Z^ 
< mV cosy > 


. . V Y 

sin <p cosy cos x tan— 


+2<a^tan y sin % cos— 
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h = V siny 


+ VRcosycosx 

X ~ Y 

rcos — 

R 

^ _ VRcosy sinj 

r 

The following are the energy-state approximations of these equations. 
( 1 )' 

m = -CT 

X = V 

mg 

L- mg 

(2)' 

m - -CT 
X = V cos% 

Y = V sin x 

E = P 
L sin 6 

x = — — 

mV 

Lcos<j> = mg 
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(3)' (with T = 0 and m = const ) 


E = - 


VD 

tog. 


L e V 

O = cos^ — — + — +2<u cos ^ cos 

zn V V r 


L V y . Y 

X~ sin^ cos ^ tan 2a> sin — 

mV r R R 


y VRcosx 

X Y 

r cos— 

R 


Y = 


VRsinz 

r 
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